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Magnetic properties of Mn-doped GaAs are re-investigated witiiin a realistic multiband description 
of the host valence bands. We explicitely demonstrate that the recent Monte Carlo (MC) simulations 
performed on a large scale supercomputer suffer from severe shortcomings. Indeed, it is shown, using 
identical parameters that (i) the calculated Zeeman-splitting largely underestimates that obtained 
from first principle studies, (ii) the couplings exhibit strong RKKY oscillations, (iii) the stability 
region for ferromagnetism is much narrower than obtained previously and (iv) the calculated Curie 
temperatures appear to be at least one order of magnitudes smaller. We show that the proposed 
choice of physical parameters cannot describe the physics in (Ga,Mn)As. 
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In the recent years, dilute magnetic semiconductors 
(DMS) have risen a lot of interest because of their pos- 
sible use for spintronic devices [l| (e.g. a spin-polarized 
light-emitting diode, etc.). A key goal is the achievement 
of a Curie temperature Tc well above room tempera- 
ture and a proper theoretical description of the involved 
physics was mainly controversial. From the theoretical 
point of view there are mainly two different kinds of ap- 
proaches, (i) realistic bandstructure model studies (k-p, 
Kohn-Luttinger, empirical tight-binding) ^ and (ii) so- 
phisticated material specific ab-initio based (LSDA, SIC- 
LDA, ...) calculations [3]. Note that the prediction of a 
large Tc ~ 700 K in GaMnN !4], has triggered a con- 
siderable amount of both experimental and theoretical 
work. Later it has been demonstrated that this huge Tc 
results from several drastic approximations namely (1) 
perturbative treatment, (2) mean-field (MF)and (3) vir- 
tual crystal approximation for disorder [3, |6| . Recently 
it was shown that the essential properties of DMS can 
already be described within a single-band model [7|, [Sj . 
Therefore a model with a realistic bandstructure and a 
non-perturbative coupling between carriers and magnetic 
impurities should provide a powerful tool to study both 
transport and magnetic properties in the wide family of 
DMS. 

Recently such a study was performed for the case of 
GaMnAs using large scale Monte-Carlo (MC) simulations 
[Ol-KLlj. Although a MC approach can be considered as 
essentially exact, the study of dilute magnetic systems 
of classical spins coupled to quantum carriers should be 
done in a very careful way. Our intent is to discuss and 
compare our calculations based on a two step-approach 
(TSA) to the MC simulations of Ref.[i3l. As it will be 
shown several crucial features are incorrect in Ref.[10| 
and explanations will be provided. 



We start from the following V-J model Hamiltonian 



H = 



V t'^^c^- c 



'jl3a 



Jpd^PiSiSi + V^ 



Pirii. 



ij^a/S^a 



(1) 

The first term provides a realistic multiband description 
of the host valence bands. Here c|^^ (ciacr) is the cre- 
ation (annihilation) operator of a carrier with spin a at 
lattice site i. a, (3 denote the three different p-orbitals. 
The hopping matrix elements i"- are restricted to nearest 
neighbor sites and are determined from known Luttinger 
parameters 7i [l^l. We used 71,2,3 = {6.85,2.1,2.9} as 
used in Ref.[lO| for GaAs. The second term corresponds 
to the interaction between itinerant carrier spins s^ and 
localized impurity spins Si (S = 5/2 for Mn^+). Jp^; is 
the local p-d coupling. The last term (on-site impurity 
potential scattering term) will be neglected (V = 0) as 
done in Ref. 10]. The relevance of the on-site term will 
be discussed in what follows. The parameter pi takes the 
values Pi ^ I ii the site i is an impurity site otherwise 
Pi = 0. Now the general procedure is as follows: First 
we calculate the total spectrum of the Hamiltonian de- 
fined in Eq. ([T]) for a given set of disorder configurations 
assuming the Mn^+ impurity spins to be fully polarized. 
Then the Mn-Mn effective exchange couplings Jy (see 
Ref. ISl, |lj| for details) between two impurities located 



at Ti and r, are obtained by the relation: 



Jij — 



AttS^ 



Jm 



f{iu)Tr^{j:M{u^)tM{u^)}du. (2) 



Note that Eq.(I2]) depends on the impurity concentration 
X and on the hole density p. In our model the self-energy 
reduces to E^ = Jp^S-l and Gfj{uj) is the Green's func- 
tion of the system while the trace is taken with respect to 



the orbitals. f{uj) 



_ (pI3(u-Ef) 



1 



is the Fermi func- 



tion and Ei? is the Fermi energy. Finally the disordered 
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FIG. 1. Empirical tight-binding (ETBM) density of states 
for Gai_2,Mna;As for x = 0.085 and different values of JpdS. 
The dashed lines indicate the chemical potentials for T — 
K corresponding to the hole density p — O.TSx. 
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FIG. 2. Zeeman splitting AEv{x) as a function of the Mn con- 
centration X for Gai_2,Mna;As calculated within LSDA [i][la] 
[ii],17.] and within the V-J model. 



Heisenberg model, 



H 






(3) 



defined with the effective exchange couplings of Eq.Q 
is solved within the self-consistent local random phase 
approximation (SC-LRPA). This method [1, 0, [11 has 
been proven to be a reliable tool to calculate e.g. Curie 
temperatures, the magnon excitation spectrum and op- 
tical conductivity. The TSA was able to reproduce both 
ab-initio results and experimental data. 



RESULTS 

For a fixed Mn impurity concentration of a; = 0.085 
in Gai_a;Mna;As the density of states (DOS) is shown in 
FIG. [TJ for different values of the effective p-d coupling 
strength Jp^S. These results are in perfect agreement 
with Ref.[lO|. The DOS is almost unchanged compared 
to that of the pure system for Jp^S < 3 eV. This already 
gives a first indication/hint that this range of parameters 
corresponds to the perturbative RKKY regime, this will 
be confirmed in what follows. Let us proceed further by 
calculating the Zeeman splitting /S.E^ (x) = E}^^^ — E^^^ 
i-^max is ^he largest eigenvalue in the corresponding a 
sector) as a function of the Mn concentration. It is de- 
picted in FIG. [2] together with available ab-initio results. 
For JpdS = 1.2 eV A£'„ is found to be in perfect agree- 
ment with the mean-field expression AE^^^{x) = xJpdS, 
the accordance is still reasonable for Jp^S = 3 eV. Re- 
mark that, it is now wi dely accepted that Jp^ « 1.2 eV 
in Mn-doped GaAs la, llflf . Thus the realistic value of 
JpdS should be 3 eV since S = 5/2 for Mn^+. However in 
Ref.[10| JpdS = 1.2 eV was used, which is almost three 
times smaller than the realistic value 20|,|21|. From FIG. 



[21 one can clearly see that for both Jp^S = {1-2, 3} eV the 
calculated Zeeman-splitting largely underestimates that 
obtained from first-principle studies [16|, |l7[ • For exam- 
ple, for X = 0.05, one finds AEy = 0.06 and 0.2 eV for 
Jpd — 1.2 eV and 3 eV respectively, in contrast to 0.65 
eV obtained from LSDA calculations. Thus, even the 
correct value of JpdS does not lead to an agreement with 
first-principle results. In addition, experimental studies 
22l426| and ab-initio LSDA calculations tl7| indicate the 



existence of an acceptor level or bound hybridized Mn pd- 
state at Ej, « 112.4 meV above the valence band. This 
impurity state is absent for both values of JpdS = {1.2, 3} 
eV in the limit of a; — >■ 0. In order to recover the correct 
Zeeman-splitting and impurity acceptor level energy us- 
ing JpdS = 3 eV, one has to include a finite additional 
impurity potential scattering term V (see Eq.(IT])) [2'/|. 
Within a single-band model this term has been shown 
to be a crucial ingredient to understand magnetism and 
transport properties in III-V compounds [71, Isj. Hence, 
our conclusion contradicts that of Ref . [9| , that a finite V 
is irrelevant. Let us now discuss the nature of the mag- 
netic couplings obtained within our multiband model. 
The Mn-Mn exchange couplings Jij are calculated ac- 
cording to Eq.(15]). In comparison to the MC studies of 
Ref. lOj our calculations are carried out for much larger 
systems (2048 lattice sites vs. 256) and the average over 
the disorder is performed over up to 800 vs. about 5 con- 
figurations in the MC simulations. Thus we expect our 
results to have much less finite size effects and reliable 
statistics. For different values of JpdS, J{R) ■ R^ a,s a 
function of R is shown in FIG. [31, where J{R) = {Jij)dis 
and R = jr^ — r^j. Our results show clearly a major differ- 
ence in the range, magnitude and nature of the effective 
exchange couplings. For small values of JpdS={1.2, 3} 
eV J{R) exhibits undamped long range RKKY oscilla- 
tions. Hence, the couplings obtained for these values 
of Jpd are inconsistent with those calculated from first 
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FIG. 3. Mean effective exchange couplings J{R) ■ R^ as a 
function of distance R {a is the lattice constant) for several 
values of Jp^S and p = 0.75a; using x = 0.085. Please note the 
oscillating RKKY nature of the couplings for JpdS = {1.2, 3} 
eV in contrast to the short-ranFR ferroma.p'netic counlinp's for 
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FIG. 4. Mean effective exchange couplings J{R) ■ R? /{J-i ■ -Rf) 



as a function of distance R for several values of Jp 
p=0.75x using x — 0.085 in the [110] direction. 
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principle studies for which RKKY oscillations are ab- 
sent (see Ref. [3[ ) . This disagreement indicates that either 
the choice of parameters used to describe the physics of 
GaMuAs is incorrect or the model itself is inappropri- 
ate. In contrast, in the strong coupling regime Jp^S = 
{8, 16} eV, the exchange integrals appear to be of short- 
range ferromagnetic type. This regime can be seen as 
a precursor of the double-exchange regime, JpdS = cxd. 
Both the RKKY and exponentially damped nature of 
the exchange couplings can be clearly seen in FIG. ID 
Note that for Jp^S < 3 eV the Hamiltonian actually de- 
scribes the physics of II- VI Mn-doped systems as, e.g. 
CdMnTe or ZnMnSc. Thus, for Jj,dS < 3 eV long-range 
ferromagnetic order is unlikely |28| , because of high frus- 
tration. Instead one would expect a spin glass phase. 



Finally, in FIG. [5] we present the calculated Curie tem- 
peratures Tc obtained by diagonalizing the Heisenberg 
Hamiltonian Eq.ljS]) within the SC-LRPA using different 
hole densities. For comparison both T^^^ (from Ref.[10|) 
and the mean-field virtual-crystal value (MF-VCA) T^^^ 
are shown. If we first focus on the largest hole density, p 
= 0.75x, we observe that T^'-'' is surprisingly larger than 
rpM_F ^ ^jj^jg jg jjj^ contrast to the expectation that the MF 
value should overestimate the real critical temperature. 
Even in the region of small JpdS (< 3eV), the MF value 
is zero whilst the MC value is always finite and relatively 
large. E.g. for Jp^S = 3 eV the MC value is about 500 K. 
Therefore one can question the reliability of the present 
MC calculations. On the other hand, the SC-LRPA Tc 
is always much smaller than the MF values for the whole 
parameter range of Jp^S. Ferromagnetism is possible only 
for sufficiently large values of Jp^S. We find that the crit- 
ical value is about 5 eV for p = 0.75x and about 2 eV 
for p — 0.3a;. Our calculated Tc are approximately one 
order of magnitude smaller than those obtained from MC 
simulations. For example for Jp^S =12 eV we have ob- 
tained 400 K within SC-LRPA whilst the MC simulations 
value is 2300 K. Note that for lower carrier density the 
critical temperatures are smaller but the ferromagnetism 
already appears for smaller values of Jp^S. This is ex- 
pected since the effect of disorder is stronger at lower 
carrier concentration and RKKY oscillations will be sup- 
pressed for smaller values of JpdS. Additionally, at lower 
density the period of the RKKY oscillations is larger thus 
the frustration effects weaken when carrier concentration 
is reduced. Let us now explain the origin of the disagree- 
ment between our calculations and the MC simulations. 
The system considered in the MC study consists of only 
256 lattice sites with typically only 20 localized spins 
and about 15 carriers/holes in the whole cluster. In ad- 
dition the statistical average was carried out over about 
five configurations of disorder only. As it was already 
pointed out in Ref. 29], the Tc is strongly size depen- 
dent and huge fiuctuations of the critical temperatures 
distribution were observed for such small system sizes. 
For JpdS < 4 eV, the smallness of the cluster used in the 
MC study does not resolve the asymptotic RKKY tail for 
JpdS = 1.2 eV leading to finite and large Curie tempera- 
tures. In the large coupling regime, Jp^S > 8 eV the T^^'^ 
overestimates the real critical temperatures both due to 
insufficient statistical sampling and finite size effects. On 
the other hand, it is clear, that these essential numerical 
requirements are difficult to fulfill within standard MC 
calculations. At last, one should underline, crucial differ- 
ences between the two-step approach and MC simulations 
concerning the way Tc is determined. Within SC-LRPA 
Tc is directly calculated from a semi-analytical expres- 
sion whilst in the present MC study it is extracted from 
the temperature variation of the magnetization curve. A 
more accurate way should be to use Binder cumulant in 
order to avoid additional errors, for the cost of signifi- 
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FIG. 5. Calculated Curie temperature To for Gai_a;Mna:As 
with X — 0.085 using the V-J model in connection with V 
= 0. The SC-LRPA calculations were done with a system of 
4- (22)"^ sites (3600 impurities) and the average was performed 
over 100 disorder configurations. 



cant computing time and memory needed for the finite 
size analysis. On the other hand, the SC-LRPA, allows 
the use of very large systems, typically of the order of 
4 • (30)'^ thus containing about 10"* impurities. 



CONCLUSION 

In this work we have studied the magnetic properties of 
Gai_2;Mna;As including a realistic bandstructure of the 
host material. We have demonstrated that the param- 
eters used in Ref . [10| . namely Jp^S = 1.2 eV and V = 
0, cannot describe the magnetic properties in GaMnAs. 
The appearent agreement pointed out by those authors 
results from several levels of "approximations". Indeed, 
the direct comparison between our two-step approach 
and the recent MC simulations revealed several short- 
comings. The MC study suffers from a) finite size effects, 
b) an insufficient statistical sampling and c) a less accu- 
rate and approximate procedure for the determination 
of the Curie temperature. It has been shown, that the 
critical temperatures Tq'^ are largely overestimated for 
these reasons and even larger than mean- field values. We 
remark, that though the MC approach is essentially ex- 
act, the study of dilute magnetic systems has numerical 
requirements hard to fulfill by standard MC techniques. 
Furthermore we point out the necessity of including a fi- 
nite V (impurity potential scattering term) in the multi- 
band V-J model in order to describe the ferromagnetism 
in (Ga,Mn)As properly. 
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